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Abstract 

^ ■ We present an extension of the Standard Model (SM) based on the discrete flavor 

(N 

00 



symmetry ^3 which gives a neutrino mass matrix with two-zero texture of Fritzsch-type 

■^ ■ and nearly diagonal charged lepton mass matrix. The model is compatible with the 

• ■ normal hierarchy only and predicts sin ^13 ~ 0.01 at the best fit values of solar and 

C^ • atmospheric parameters and maximal leptonic CP violation. 

O' 
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1 Introduction 

Although there is a robust evidence that neutrinos are mixed, many aspects of the neutrino 
physics are not clearly understood yet. Among them, the comprehension of the values of 
the masses and mixing and the differences with respect to the quark sector are an open 
problem whose solution seems to be quite far from being found. Recent data from neutrino 
oscillations produced the following results: 

0.36 < sin^ 023 < 0.67 0.27 < sin^ 9^2 < 0.38 sin^ ^13 < 0.053 , (1) 
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and 

9 n? V in-3 pT/2 < Arr)^. < 9 7R v in-3 pT/2 

(2) 



2.07 X 10"=^ eV^ < Amlt^ < 2.75 x lO'^ eV^, 



7.03 X 10-5 eV^ < Am'l^i < 8.27 x lO'^ eV^, 

at 99.73% confidence level [Ij (see [2] for other recent interpretations of the neutrino data). 

We have only hints coming from cosmological observations that the absolute values of 
the neutrino masses should be less than 1 eV ^. In the quark sector the situation is quite 
different: not only the masses and the hierarchy in the up and down sectors are better 
known but also the mixing angles are well measured and strongly differ from the neutrino 
ones. A successful ansatz to reproduce these features in the quark sector is the Fritzsch-like 
texture |4], where both the up and down quark mass matrices have a simple form 

/ A \ 



M= A* C B 
\ B* D J 

Such a matrix (already described in, e.g., [5]) gives the well know relation 



(3) 



tan^i2 = W— , (4) 

V m2 

which predicts the Cabibbo angle whose small value is a consequence of the strong hierarchy 
in the masses. A texture as in eq.Q can also be employed for the Majorana neutrino mass 
matrix; this is a particular case of the class of two-zero texture [0] which, together with the 
two relations Amatm = ""^i — f^i and Arusoi = Tn\ — m\ , fix the absolute neutrino mass scale 
as suggested in [7]. Unlike the quark sector, the solar and atmospheric angles can be large 
due to the fact that in the neutrino sector the hierarchy is not so strong. 

Although a vast class of Fritzsch-like textures (and their phenomenological consequences) 
has been already studied in the literature, in this paper we propose a leptonic model based on 
the permutation symmetry S-^ which naturally gives rise to a Fritzsch-type neutrino Majorana 
mass matrix (and, in addition, to a nearly diagonal charged leptons). At tree level, the tau 
lepton acquires a mass via the spontaneous electroweak symmetry breaking (ESB) driven 
by one 5*3 doublet and two S^ singlets, whereas the electron and the muon remain massless. 
Higher order operators, mediated by just one Standard Model (and ^3) scalar singlet (called 
the flavon) are responsible for me, m^ 7^ 0. In the neutrino sector, the Majorana mass matrix 
is generated by dimension five |8J and six operators. 

The paper is organized as follows: in the next section we introduce the model; the scalar 
potential is studied in SecJS) the lepton and neutrino mass matrices are introduced in SecsJH 



and El respectively whereas their phenomenological consequences are discussed in SecEl 
Seed is devoted to our conclusions. 



2 The model 

We propose a model based on S3, the group of permutations of three objects, which is the 
smallest non-Abelian discrete group. S3 contains one doublet irreducible representation and 
two singlets. This feature is useful to separate the third family of fermions from the other 
two and has been already used for model building [9] . For pioneers papers see [10] (and also 
references in [TT]). 

The group 5*3 has two generators S and T satisfying the following relations: 

S^ = T^ = (5T)2 = 1 . (5) 

One possible realization is the so-called "T-diagonal" basis where 



s= T= z :, . (6) 





with Lo' = e*^/^'^. The tensor products involving pseudo-singlets are given by 1' x 1' = 1 and 
1' X 2 = 2 while the product of two doublets is 2x2 = 2 + 1 + 1' which, in terms of the 
components of the two doublets A = (ai, 02)""" and B = {pi, 62)""" in the T-diagonal basis, are 
as follows: 

0162 + 0261 G 1 / 0262 ] ^ 2 m 

aib2 — a2&i G 1' I aibi J 

The product 2* x 2 is similar to 2 x 2 with the exchange of ai f-)- 02. 

Construction of the model 

The Higgs sector is extended from one S'f/(2)2,-doublet to two S'f/(2)L-doublets, Ho = 
{Hi, H2) belonging to a doublet irreducible representation of 5*3 and other two SU{2)l dou- 
blets, Hs and Hg, belonging to singlet representations of 5*3. We also introduce an elec- 
troweak scalar singlet x which turns out to be relevant to give a non-vanishing electron 
and muon masses. In order to have nearly diagonal charged lepton mass matrix we assume 
two further parity symmetries, so that the global discrete symmetry group of the model is 
G = S30 Z^^ Z2. The matter assignment under G ® SM is summarized in TabiH 
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Table 1: Matter assignment of the model, to is the Z^ charge u = e*^/^'^ and Y is the SM 
hypercharge in the convention Y = 2{Q — T3), where T3 is the third component of the SM 
SU{2) doublets. 

3 The scalar potential 

The most general Higgs potential invariant under G x SM is as follows: 

+ {\2HlHn + \3HlHs + KH'^H'sM' + \,[{HlHn)? + Ae [(if l,//z>) if 

+\',{HlHnhH'^H's + \'i{H^^H'sUH'^Hn), + K'{{H^nH's)l + /^•c.)+ (8) 

+K,{HlHn)iHlHs + \UH^j,H'sUH'^sHd)2 + A^ ((i^l^iis)! + /i.c.) + 



rt H-t 



+\niHj,Hj,UHDHD)2 + \i2{H'^H'sf+ 

+Xi3'H'^H'sHIHs + Xn"{H'^H'^HsHs + h.c.) + Xn"'H'^HsHlH's 

where we used the subscripts 1, 1' and 2 to refer to the S3 contractions when necessary and, 
for any Higgs fields, H = —ir^H*. In the case of real vev's, that is 



{Hd) = {vi, V2), 
(Hs) = vs, 
(H's) = v's, 



(9) 



the potential can be written ai 

V = ( An + As + Ae) K + 4) + 4^8 + v'^Xu + v^v'^Xu + v'^fn + vlfi3+ 

+ {fi4 + vlX3 + v'iXi)x'' + Xix^ + {vlXio + v§Xg + fi2 + X2x'){vl + vf)+ (10) 
+ {2vy^{X5 - Xe + X7)) . 

The minima of V are found solving the minimizing equations: 

1^ = 2t;i [X2VI + Xiovl + Xgv'l + 2t;2(An + A5 + Ag) + 2vl{X^ - Ag + X^) + /is] = 0, 

g = 2^72 [A2t;^ + XiQvl + XWs + 2t;i(Aii + A5 + Ag) + 2t;2(A5 - Ag + A7) + /is] = 0, 
fj = 2vs [X^vl + Aio (t;? + vl) + 2t;|A8 + v'lXi, + /is] = 0, 
g = 2< [A4t;2 + A9 {vl + t;i) + 2t;f A12 + t;iAi3 + /ii] = 0, 

fl = 2^;^ [/i4 + 2Ait;2 + A2 {vl + t;i) + Agt;! + X^v'l] = . 

(11) 
The second equation is satisfied for f 2 = 0. From the remaining equations we can easily get 

the vevs of the other scalars in terms of the couplings of the Higgs potential; in particular, 

a solution with i^i 7^ can be found and the vev alignment of the 6*3 Higgs doublet assumes 

the structure: 

{Hn) = {v,^). (12) 

For this vev configuration, it is possible to find a huge region of the Higgs parameter space 
where the eigenvalues of the Hessian of the potential are all positive and therefore where 
the Higgs potential has a local minimum. Note that a solution of the form {Hd) = (0, f ) is 
physically equivalent to eq. (iT2l) . producing the same phenomenology in the charged lepton 
and neutrino sectors. In fact it corresponds to the exchange of Li with L2. We also verified 
numerically that, in the large parameter space where eq.f lT2|) is a minimum, other solutions 
like (Hd) = V (1, 1) do not produce positive definite Hessian. The mass spectra of the Higgs 
particles will be discussed elsewhere. 



*Where Ai = A^ + A^' + A'/' for i = 9, 10, 13. 



4 Leptons 



The most general Lagrangian invariant under G x SM is given by: 



-j-LdHdIrdX* + -j-LdHsIroX + VsL^HsIr^ 



(13) 



where A is the cut-off scale. Higher order terms only appear at 0{1/A^) and will be consid- 
ered negligible for our discussion. From eq. flT2|) the charged lepton mass matrix is: 



Ml 



( f^sv^ 








>^x 


>sv^ 





\ 





V'iVs 



(14) 



When v^ is equal to zero only the r lepton is massive. The electron and muon masses are 
generated by the vev of the scalar x ^i-nd are then suppressed by the large scale A. The 
matrix MiMJ has three distinct eigenvalues that can be identified with the squared charged 
fermion masses as: 

r-2 



m„ 



.2„.2 



.2„.2 



m 



m. 



v'y'i + 2vlyl - vyi^Jv^yl + Avlyl 
y Wvl + '^vlyl + vyi^Jv^yl + Avlyl j 



(15) 



2 2 



where we introduced the short-hand notation e = v^/A. We see that for e ^ 1, the hierarchy 
among the r and the lightest charged leptons is easily reproduced although the latter, in 
absence of any fine-tuning among the Yukawas and/or the Higgs vevs, are expected to be 
of the same order of magnitude. We address this question in the next section. The mass 
matrix for the charged leptons can be written in terms of the physical lepton masses as: 
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y/rriem^ 
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and the squared matrix MiM^ is then diagonalized by: 



m. 



(16) 
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(17) 



5 Neutrino 

The neutrino masses are generated by non-renormalizable operators of dimension 5 and 6 
invariant under the group G x 5'iWlj: 



A-C = yl{LDLD)^{HDHD)l+y2{LDLD)2{HDHD)2+yl{LDHD)l{LDHD)l + 

y'i(LDHD)v{LDHD)v + yl{LnHD)2{LDHD)2 + 

yl'(LDLD)i{HsHs)i X* / ^ + y',"(LDLD)i{H'sH's)i X* / ^ + (18) 

y^L^L^{HDHD)i X/A + yl'L,L,HsHs + yl" L.L^H'sH's + 

yl{LDHD)iL2,H'g, 

where we assumed that the large energy scale which suppresses these operators is of the 
same order of the cutoff scale A. The only operators of dimension six are those proportional 
to ?/g = y^' + y^" and y^. From eq.( IT9|) the neutrino mass matrix is as follows: 



M, 



/ 2y^,{v1 + v'i)vJA \ 

2 y'M + v§)vJA {y^, + yi; + ylW y^vv's 

\ y^vv's yl{vl + v'^) ) 



/ b 0\ 

b a c 
\0 c d J 



(19) 



where y^ = y^' + yg". Before discussing the phenomenological consequences of such a matrix, 
it is useful to get an estimate of the relevant Yukawa parameters and a relation among the 
vevs V and vs- Comparing eqs.l TT^ and ( TT6l) and using the parameterization in eq.f llQp we 
get: 



m^\ 1 - 



m 



m 



a 
b 



y2Vse 



yive 
VsVs 



where we assumed that y2 + ys + y^ = y'^- We assume 



f ^ f « ~ 



O(IOO) GeV, 



e ~ 



0(10- 



^Dimension 7 operators can be built, for instance, adding the singlet x^ or doublets H'' H to the previous 
d — 5 operators and will then be neglected. 



7 



n 



(20) 



then if 

1/1 ~ 0(10-^), 1/2 ~ C(10-2), 1/2 ~ 0(10-2) 

we have the correct charged lepton mass hyerarchies. We then consider the b/a ratio 

a y" f 2 vs 

We numerically verified that (|) ~ 0{1) so the Yukawa parameters must satisfy 

With these assumptions, the hierarchy in the charged leptons is recovered, higher order terms 
with more that one fiavon insertions can be safely neglected and the largest vev is generated 
by the Ss singlet Higgs Hs that can be identified with the Standard Model Higgs. 

The mass matrix in eq. (fT9l) depends on five real parameters, one of which is related to 
the Dirac phase. The other four parameters can be fixed using the experimental information 
from both solar and atmospheric sectors, namely the solar and atmospheric mixing angles 
and squared mass differences. The model allows for correlations among the angle 615 and the 
CP phase 6 that can be easily obtained using the zeros of the Fritzsch texture. The previous 
mass matrix is diagonalized by a unitary matrix W^ as 

f/'^^M,f/'' = diag(/xi,/i2,/X3) (21) 

where /ij = rriie^^^ and 0j are Majorana phases. Writing V^ in the CKM-like forr 






U" 



i&v rM r,V „U „V „i(5, 



V „V r,V 



23-^12 '-12'^ ''13''23 '-12'-23 '^ ''12''13''23 "-13''23 



(22) 



\ ^12^23^^ ■'IS "T ■^12''23 ''23'^ ■^12'^13 ''12''23 ''13''23 / 

and using the fact that the elements {My)\\ and {My)x-i are zero (see eq. (TT9!) ). eq. (!22|) implies: 



cos ej^j (- cot ei^ cos e^3+sin6»5;'2 sin6»5;'3 sin 6*^3 e'"*'^ ) 
/^2 = A'-l cos 6*^3 sin 9^2+cos d\^ sin d\^ sin 6^3 e'*-- ' 

_ _ cos e\^ cos^ e\^ sin 61^36-'^^ 

^^ ~ ^-'^sin6IJ'3(cosej3sin6»J'2+cose5'2sine5'3sin9^3e»''i') 



(23) 



Our model is compatible with the normal mass ordering only because the ratio |/X2p/|/^3p is 
always less than 1; expanding it up to second order in sin6'5'3 we get: 

'^''' - cot^ Q\^ cot^ ^-3 sin^ Q\^ + 0{st) , (24) 



l/^sP 



^ We have used the short-hand notation sf, — sin 0f,- and cj', — cos flf, . 

*j ^j ''J ^j 



and we checked that higher order corrections do not modify our statement. The mass 
differences are written as: 



Aml^i = ml 



cos ^23 (cos 0J'2 cos 9^^ csc^ 9'(^ - 2 cos 6" cot 9'(^ sin 9'(^ sin 9'^^) 



(25) 



and 

A 2 2f. cos^ ^1^2 cos^ ^5^3 cot^ 9'(. sin^ 9^. \ 

Aml^ = ml{^-l + ^ ^- ^ 23 j ^ (26) 

where 

D^ = cos^ ^^3 sin^ 9'(^ + sin 2^;;'2 sin 2^^3 sin 9'^^^ cos ^72 + cos^ 9'^^ sin^ ^^^3 sin^ ^^3 . (27) 



From the ratio a = Aml^J Am\^^ we find a relation between a and the mixing angles s 



12' 



-Sis and S23 



a = 2 sin^ ^J'3 cot ^^3 esc 9^^ sec ^^2 (- sin 9^^ cos 5^ + cot 2^J'2 cot ^23) • (28) 

which will be used below to constrain the physical 9i^ and 5. 

6 Phenomenology 

To study the phenomenological implication of our model, it is necessary to relate the param- 
eters in eq. (l22p to the physical ones. This can be achieved introducing the rotations from the 
charged lepton sector described in SecJH the resulting mapping is a set of implicit relations 
that are quite cumbersome and will not be explicitly presented here. We limit ourselves to 
describe the procedure which allows us to extract the predictions of our model. The lepton 
mixing is defined by V^ = UjU'^ and we can write 

U^ = UlV, (29) 

where V is parametrized in the standard form as in eq. 0221) replacing s\, and cjl- with the 
physical (that is measurable) Sij = sin 9ij and Cij = cos 9ij, respectively. Taking the ratio of 
U23 and [733 from eq. (12 9 p we find an expression for S23 in terms of the physical parameters 
9i2, 923, ^13 and the Dirac phase 6 (and the corrections from the charged leptons). In the 
same way, always using eq. fl29]) . we can express s^j'g and 5^3 as a function of 9i2, 6*23, ^13, 6; 
finally, 6^ is the argument of the element (13) of the matrix UlV. In this way we have all 
the parameters 9^2, ^13, ^23 ^^cl 61, as a function of the neutrino mixing angles 6^13, 6*12, 6*23 

9 



and the phase S. These relations can be inserted into eq. (12 8 p to get an imphcit connection 
among the mixing parameters and a, which is a characteristic of our model. Also the lightest 
mass eigenstate can be related to the same parameters and Am^^,^ using eq. ( I26l) . 

In the left panel of Figil] we show the dependence of sin^ ^13 as a function of 6 taking 
^12, ^23 and a inside their experimental ranges. In particular, the solid line represents the la 
correlation when also the other parameters are left free to vary in their la allowed ranges 
quoted in [1], whereas the 2a correlation is represented by the dot-dashed line. Finally, 
dashed line is the relation obtained when 6'i2,6'23 and a are fixed to their best fit values. 
We also included the upper limit on sin^ 613 at 3a (upper horizontal dashed line) and the 
best fit value of ref. [T| (lower horizontal dashed line). We can see that, even considering 
the 2a uncertainty, the predicted values for sin^ ^13 are different from zero so that, to a very 
good accuracy, our model is compatible with deviation from ^13 = for any value of the 
CP violating phase. The precise value of 6*13, however, relies on the assumed magnitude for 
6; in particular, the CP conserving case 5 = is the most promising one to allow large 6*13 
(even above the current limits) whereas around 6 ~ ±7r we get the smaller ^13 allowed in our 
model. It is interesting to observe that, in the case of maximal CP violation \j and for the 
other oscillation parameters to their best fit values, the predicted sin^ ^13 is fully compatible 
with the best fit value obtained in pTj, sin^ ^13 ~ 0.01. Notice that, in the case of diagonal 
charged lepton mass matrix, the pattern of the 613 — 6 correlation would have been quite 
similar, as it can be seen investigating the right panel of Fig{Tl The fact that the corrections 
coming from Ul in eq. flT7|) are as large as the values of sin^ 613 is responsible for lowering 
the allowed Oi^ for 6 ~ ivr. For maximal CP violation at the best fit point, the Jarlskog 
invariant [Hj is as follows: 

J = C12S23C13S12S23S13 sin 6 = 0.023 . (30) 

The next observable we want to discuss is the effective mass rriee entering in the neutrino- 
less double beta decay. In the basis where the charged leptons are diagonal, rriee is nothing 
but the (11) element of the neutrino mass matrix. According to eq. flTQl) . this should vanishes 
as long as the rotation in the charged leptons is proportional to the identity matrix. Since 
this is not the case, a non- vanishing mee is generated by the rotation ( 1T71) and it is expected 
to be small because of the smallness of its off-diagonal entries. This is what we can observe 
in Fig. (12]), where we plot the model predictions for TTige as a function of the lightest neutrino 
mass nil. For rrii below (9(10~^) eV we get |T7iee| ~ 10^^ eV and then outside the range of 



''Maximal CP violation can be observed in incoming experiments T2K and NO^^A, see for instance [T^13) . 

10 



future experimental sensitivities. We also see that the allowed range for the lightest neutrino 
mass is around 10^^ — 10~^ eV; this is because, as already mentioned in the introduction, 
the Fritzsch texture gives a correlation between 9^2 u and the ratio mi/m2 that, together 
with the two measured square mass differences, fix the absolute neutrino scale in this range. 





Figure 1: Left panel: correlation among 6 and sin^ ^13 as obtained in our model. The la 
result, obtained varying the other oscillation parameters also in their la allowed ranges, is 
showed with solid lines, whereas the 2a result is showed with the dot-dashed line. The dashed 
line is the relation obtained when 612, 623 and a are fixed to their best fit values. Horizontal 
lines represent the upper limits on sin^ ^13 (upper dashed line) and the best fit values (lower 
dashed line) from |^. Right panel: the same as the left panel but assuming exactly diagonal 
charged lepton mass matrix. 



7 Conclusion 



In this paper we have studied a leptonic model based on the discrete 5*3 permutation flavor 
symmetry. We extended the scalar sector of the Standard Model by introducing three more 
Higgs doublets and one scalar singlet. We have carefully studied the problem of the mini- 
mization of the potential and found that a solution of the form (f , 0) for the Higsses in the 



^Note that in the case of diagonal charged leptons the angle 6*12 corresponds exactly to 9^2 and therefore 
from eqs. ^ and (P^. the solar mixing angle does not depend on the absolute scale of neutrino mass mi, 



while in our case this relation acquires a small correction proportional to \Jme/m^ 
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Figure 2: Model predictions for |mee| os a function of the lightest neutrino mass mi. We also 
show the allowed regions for the normal hierarchy (gray hand). The two dashed horizontal 
lines represent the experimental sensitivity of some of the forthcoming experiments while 
the dashed vertical line is the upper limit for the sum of the absolute neutrino masses from 
cosmological data. For references to experiments see /75l - F7^ . 

5*3 doublet representation is a viable minimum of the potential. With such a minimum, we 
obtain a two-zero Fritzsch-texture for the neutrino mass matrix and a nearly diagonal and 
hierarchical charged lepton mass matrix. As a consequence of the two zeros of the Fritzsch 
texture, we get a strong correlation between the reactor angle 9i^ and the Dirac CP phase 
5. In particular, for 5 ~ ±7r/2 we predict sin^ ^13 ^ 0.01, a value which is very close to the 
best fit value quoted in [1]. Beside the reactor angle, we also investigated the prediction for 
the effective mass m^e governing the rate of the Qvfifi decay, founding niee ~ 10^^ eV , one 
order of magnitude less than the sensitivities of the future experiments. 
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